Single ion anisotropy of Mn doped GaAs measured by EPR. 
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Electron paramagnetic resonance (EPR) study of MBE grown Mn doped GaAs is presented. The 
resolved fine structure allows us to evaluate the crystal field parameters of the spin Hamiltonian. 
The obtained cubic constant is o = —14.1 • 10"'' cm~^. The axial field parameter, D, increases with 
Mn concentration, x, i.e., with the strain of Gai_i,Mna,As layers. Extrapolation of D shows that 
the single ion anisotropy is the important contribution to the magnetic anisotropy which is observed 
in ferromagnetic layers with greater Mn concentrations. The analysis of the EPR linewidth shows 



that native defects of the concentration of 5 • 
of crystal field fiuctuations. 



10 cm , but not the Mn ions, are the main origin 



PACS Nos. 71.55.Gs; 76.30.Fc; 75.50.Pp; 75.30.Gw 



Gai_2;Mna;As is a semiconductor which exhibits ferro- 
magnetic properties . It is one of the most promising 
diluted ferromagnets which can be applied in magnetic 
semiconductor devices. The ferromagnetic phase occurs 
already for x— 0.02 and for x= 0.04 the Curie temper- 
ature reaches 90 K . Some logic structures, where the 
magnetic phase is switched on by an electric gate have 
already been built ||^. 

Ferromagnetic layers grown by low temperature molec- 
ular beam epitaxy (MBE) exhibit an axial magnetic 
anisotropy which plays a crucial role in magnetic prop- 
erties of Gai_xMna;As layers and multi layer structures. 
The sign of the anisotropy depends on the substrate. For 
a tensile strain the direction parallel to the growth di- 
rection is the easy axis while for compressive strain the 
easy direction lies in the grown plane. The magnitude 
of the anisotropy field is of the order of 5-50 mT. Def- 
initely, the anisotropy does not originate from magnetic 
dipole interactions. Such a shape anisotropy should be 
by an order of magnitude smaller and the easy axis would 
have to lie in plane. The complex nature of the va- 
lence band structure leads us to expect pseudo-dipole or 
Dzialoschynski-Moriya interactions . They could lead 
to a broadening of the resonance linewidth but their the- 
oretical estimations is not simple. To distinguish whether 
the anisotropy originates from a single ion anisotropy of 
Mn^"*" ion or from a non Heisenberg exchange coupling 
between Mn^+ ions we have undertaken a detailed elec- 
tron paramagnetic resonance (EPR) study of Mn in MBE 
grown Gai_2;Mna;As. 

Some reports on EPR of Gai_j;Mna;As were already 
published |@-^. The hyperfine structure (HFS) is com- 
monly observed but the fine structure (FS) of the EPR 
spectrum has not been resolved until now. Almclch and 
Goldstein estimated the upper limit of the cubic com- 
ponent of the crystal field parameter, a, in bulk sam- 



ples 0. The axial crystal field component, described by 
the parameter D, which may occur in strained structures 
only, has not yet been analyzed. In this paper we present 
EPR studies of weakly doped Gai_a;Mna;As compounds 
growth by MBE. We are able to find a resolved FS and 
to evaluate the single ion anisotropy and fiuctuations of 
the crystal strain. No broadening of EPR lines caused by 
non-Heisenberg component of Mn-Mn coupling has been 
found. 

The samples were grown by low temperature molecular 
beam epitaxy (LT MBE) in a KRYOVAK MBE system 
]lOt . The calibration of Ga and Mn sources allows us 
two estimate the Mn content with an accuracy better 
than 0.1%. The composition of samples with low Mn 
concentrations was verified by SIMS measurements and 
by the integrated amplitude of EPR signal. The thickness 
of low concentration samples was about 2 micrometers. 
The XRD measurements showed that all samples with 
Mn content from 10"^ to 10~^ are coherently strained 
by the GaAs(lOO) substrate. 

In EPR studies we found that samples with x < 1.5 ■ 
10"'^ are characterized by a well resolved HFS and some 
weakly resolved FS. For x > 2.5 • 10~'^ only a single broad 
line is observed. The spectra of Gai_2:Mn2;As (x = 4 ■ 
10~^) for four various directions of the applied magnetic 
field are shown in Fig. 1. The magnetic field was applied 
in the (110) crystal plane and 8 is the angle measured 
from the [001] direction. 

These spectra are characterized by g = 2 and defi- 
nitely originate from S* = 5/2 of Mn^+ ions. Schneider 
et al. found additional lines of a neutral acceptor with 
g—2.77 jl^. We observe only the signal of ionized accep- 
tors which indicates that LT MBE samples are strongly 
compensated. 
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FIG. 1. EPR spectra for various directions of external mag- 
netic field, H, applied in (110) plane. The fine structure is well 
resolved for H parallel to [001] direction (G = and T < 20 
K). The dashed line correspond to best fitted shapes. The up- 
per diagram shows the angular dependence of the EPR line 
positions. 

The six lines of the HFS structure, corresponding to 
six possible projections, Af/, of the 7 = 5/2 nuclear spin 
of Mn, are well resolved for an arbitrary direction of the 
magnetic field. The value of the HFS constant, A = 
51.6 • 10""* cm~^, agrees with that known from literature. 
The dependence of the spectrum on the magnetic field 
direction is caused by the angular dependence of the FS 
structure. The FS splitting is small and the FS structure 
is resolved only for the magnetic field directed along [001] 
(0 = 0). It manifests itself as extra lines below the first 
and above the last of the six HFS lines. 

We analyze the experimental data according to the 
standard spin Hamiltonian of the form ]|l2| : 

Hs - g^iBHS + ASI + D[S^_ - {1/3)S{S + 1)] + 
+ (a/6)[5* + Si + St- (1/5)5(3^2 ^ 35 ^ ^^j^^) 

here the first term is the Zeeman interaction, fis is the 
Bohr magneton. The second term corresponds to the 
HFS coupling which is assumed to be isotropic. The 
third and the fourth terms describe the FS structure, 
corresponding to the axial and the cubic symmetry of 
the crystal field, respectively. The parameter D is di- 
rectly related to the single ion contribution to the mag- 
netic anisotropy. The corresponding anisotropy field for 
ferromagnetically saturated spins is Ha = SD/g^s- 



The positions of 30 spectral lines, corresponding to six 
orientations of the nuclear spin, /, (M/ = —5/2, ...,5/2) 
and to five different spin flip transitions between six en- 
ergy levels labeled by Mg {Alg = — 3/2, 5/2) were 
calculated according to spin Hamiltonian (1). To eval- 
uate the spin Hamiltonian parameters we took the best 
fit of all spectra measured for various field orientations 
Q. Details of the fitting procedure is discussed below. 
The best fit was obtained for \a\ = (14.1 ± 0.3) • lO^'' 
cm~^. Within the experimental error this value does not 
depend on the Mn composition, x. The value of the axial 
component of the crystal field, D, is much smaller than 
the cubic parameter, a. For the data shown in Fig. 1, 
D = (1.3 ±0.02) • 10"'' cm-^ The dependence D on a; is 
shown in Fig. 2. 

Because the line broadening is mainly due to super 
HFS caused by the coupling of Mn^+ spins with the nu- 
clear spins of the ligands ^ , the shape of the individual 
spectral lines is assumed to be described by a Gaussian 
function which gives a much better fit as compared to 
fits when Lorentz line shapes were assumed. The satel- 
lite lines are characterized by a slightly bigger linewidth. 
In the fitting procedure we used two parameters describ- 
ing the width of individual spectral lines: the linewith 
of the central FS line, A_ffi/2,-i/2: a-nd the parameter 
describing an extra broadening of the satellite lines of 
FS: SAH = Ai75/2.3/2 - AHi/2,-i/2- The width of the 
line (3/2,1/2) is assumed to be equal to the mean value: 
Ai?3/2,i/2 = (Ai^5/2,3/2 + Ai/i/2, - 1/2)/2. The best fit- 
ted value of the linewidth is Affi/2,-1/2 = (2.0 ± 0.02) 
mT. It is found to be isotropic and independent of x. 
The additional broadening of the satellite lines, SAH, is a 
small effect. For the data shown in Fig. 1 SAH— OA mT. 
This corresponds to a fluctuation of D, SD = SAH/ gjiB: 
equal to 2-10"^ cm~^. As it is shown in Fig. 2b this value 
varies from sample to sample by about 40% but there is 
no systematic variation of SAH with Mn concentration. 

The resolved FS and HFS vanish with the increase of 
Mn concentration. To simulate the observed shape of the 
EPR spectrum it is not enough to consider the broaden- 
ing of the individual spectral lines. A good fit is obtained 
only when we treat the spectrum as a sum of the resolved 
one, originating from individual, non interacting Mn^"*" 
spins, and the unresolved contributions with averaged 
FS and HFS. With an increase of x the width of indi- 
vidual lines of the resolved spectrum (A_ffi/2,-i/2 E^nd 
SAH) practically does not change. The resolved struc- 
ture gradually disappears because of the decrease of the 
amplitude of the resolved contribution. For low x the 
width of the unresolved contribution corresponds to the 
second moment of the spectrum of non interacting Mn^+ 
spins. With a concentration increase a decrease of the 
width of the unresolved spectrum is observed. We relate 
this decrease to the exchange narrowing effect. 
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FIG. 2. (a) The dependence of the single ion axial 
anisotropy D on the Mn composition function of Mn 

composition, x. The increase of D, dD/dx — 0.09 cm~^, is 
caused by a compressive strain originating from lattice mis- 
match. The single ion anisotropy fields DS/giiB are labeled 
on the right axis, (b) The difference of the linewidth of satel- 
lite and central lines of FS, 5AH, scales fiuctuations of the 
crystal field, 5D, caused by fiuctuations of the elementary cell. 

The axial component of the crystal field anisotropy 
D is caused by a local strain acting on the ion. It 
increases with an increase of x. The data shown in 
Fig. 2 can be approximated by a linear dependence 
D = (0.09 • X + 1.05 • 10-4) cm-i. The value DS/gfiB 
corresponds to the single ion contribution to the total 
anisotropy field (when a saturated polarization of Mn 
spins is assumed). For a; = 1.5 • 10"^ the anisotropy field 
is DS/ giiB= 0.7 mT. However, since the layer strain 
increases linearly with x the single ion anisotropy ex- 
pected for x = 0.035 is DS/g^B ~ 10 mT. This value 
agrees well with the coercion field observed in ferromag- 
netic samples HJl^- It allows us to conclude that the 
single ion anisotropy is the considerable contribution to 
the total magnetic anisotropy of ferromagnetic layers. 

The fact that the extrapolation of the parameter D to 
low Mn concentration, a; ^ 0, does not tend to zero we 
relate to an additional strain which appears at low tem- 
perature caused by a difference of the thermal expansion 
coefficients of Gai-^jMn^-As, GaAs and sample holder. 

The signs of the parameters D and a cannot be di- 
rectly determined from EPR measurement. The data al- 
low only to conclude that the parameters are of opposite 
sign, a-D < 0. However, from the fact that the magnetic 
anisotropy in ferromagnetic Gai-xMn^As layers grown 



on GaAs substrate shows an easy in-plane magnetiza- 
tion, we conclude that D > 0, and, in consequence, we 
suggest that a < 0. 

Within the whole range of composition where the EPR 
spectrum is well resolved the linewidth of the central FS 
line, AHi/2, is constant. Up to x = 1.5 ■ 10^"^ we do not 
observe any additional broadening of the resolved part 
of the spectrum which could be attributed to dipole or 
pseudo-dipole (non-Heisenberg exchange coupling) inter- 
actions. The resolved part of the EPR spectrum, how- 
ever, does not originate from a statistic Mn spin but 
from those Mn ions which do not have any other Mn 
ion in the neighborhood. Because of that we are not able 
to conclude about the magnitude of a mean dipole or 
pseudo-dipole coupling between Mn ions. Moreover, the 
discussed samples are strongly compensated. Because of 
that the fact that no effect of non-Heisenberg Mn-Mn 
coupling has been found does not allow us to made any 
solid conclusion about the character of RKKY interac- 
tions. 

The additional broadening of the satellite lines of FS, 
5 AH, can be attributed to fluctuations of the local strain. 
Because the position of the central FS line is not affected 
by crystal field, the strain induced fluctuations of FS 
structure parameters, a and D (including a random di- 
rection of the local strain), do not affect the width of the 
central FS line, A i?i/2.-i/2, but they lead to a broaden- 
ing of the satellite lines only. Consequently, 6 AH scales 
the amplitude of the strain fluctuation. Assuming that 
the evaluated derivative dD /dx=0.09 cm~^ (see Fig. 2) 
is directly related to the difference of the lattice constants 
between Gai_2;Mna;As and the GaAs substrate |P,p^ one 
can estimate fluctuations of the deformation of the crys- 
tal cell, Oo, in the vicinity of Mn dopant caused by other 
impurities to be Sao = lO-'^Oo. This value is consid- 
erably bigger than the mean fluctuations amplitude as 
measured by rocking curve of XRD. The fact that Sao is 
independent of a; indicates that the fluctuations, Sao, do 
not originate from neighboring Mn ions but from other, 
non magnetic defects, probably antisites. 

Concluding, the single ion anisotropy caused by 
the layer strain is an important contribution to the 
whole magnetic anisotropy observed in ferromagnetic 
Gai_a;Mn2;As layers. Discussion of the local strain, as 
measured by broadening of the satellite lines of FS, SAH, 
allows to estimate crystal strain fluctuations. The inde- 
pendence of the crystal field fluctuations, of x indicates 
that the concentration of native defects is independent 
of the Mn content. Assuming that the range of x where 
the well resolved EPR spectrum of ionized Mn acceptors 
is observed corresponds to the range of insulating and 
fully compensated crystal, the upper limit of this range, 
X = 2.5 • 10"'^, allows to estimate the concentration of 
compensating centers to be 5 • 10"'^^ cm^'^. 
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